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Abstract. Let X be a Hamiltonian vector field defined on a symplcctic man- 
ifold (M, ui), g a nowhere vanishing smooth function defined on an open dense 
subset AI° of M. We will say that the vector field Y = gX is conformally 
Hamiltonian. We prove that when X is complete, when Y is Hamiltonian 
with respect to another symplectic form 0J2 defined on M° , and when another 
technical condition is satisfied, then there is a symplectic diffeomorphism from 
(M°,U2) onto an open subset of (M, ui), which maps each orbit to itself and 
is equivariant with respect to the flows of the vector fields Y on Af° and X 
on M. This result explains why the diffeomorphism of the phase space of the 
Kepler problem restricted to the negative (resp. positive) values of the energy 
function, onto an open subset of the cotangent bundle to a three-dimensional 
sphere (resp. two-sheeted hypcrboloid) , discovered by Gyorgyi (1968) [9], re- 
discovered by Ligon and Schaaf (1976) [15], is a symplectic diffeomorphism. 
Cushman and Duistcrmaat (1997) [5] have shown that the Gyorgyi-Ligon- 
Schaaf diffeomorphism is characterized by three very natural properties; here 
that diffeomorphism is obtained by composition of the diffeomorphism given by 
our result about conformally Hamiltonian vector fields with a (non-symplectic) 
diffeomorphism built by a variant of Moser's method [19]. Infinitesimal sym- 
metries of the Kepler problem are discussed, and it is shown that their space 
is a Lie algebroid with zero anchor map rather than a Lie algebra. 

1. Introduction. I am very glad to submit a paper for the special issue of the 
Journal of Geometric Mechanics in honour of Tudor Ratiu. I followed his scientific 
work for several years; I specially praise the superb book he wrote with Juan-Pablo 
Ortega [20]. I hope that he will find some interest in the present work. 

The Kepler problem is a completely integrable Hamiltonian system with impor- 
tant applications in the physical world: it provides a very accurate model of the 
motion of planets in the solar system, and its quantized equivalent provides a good 
model of the hydrogen atom. Moreover, several features of the Kepler problem 
make it very interesting for the mathematician: some of its motions (those in which 
there is a collision of the moving point with the attractive centre) are not defined 
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for all times, but the system can be regularized, i.e. mapped into a new Hamito- 
nian system whose motions are defined for all times; the Lie algebra of infinitesimal 
symmetries of a given energy level of the phase space depends on that energy level. 

In Section 2 some results in symplectic geometry are presented. These results 
will be used in Section 3 to explain why a diffcomorphism between the phase space 
of the Kepler problem, restricted to negative values (resp. positive values) of the 
energy, and an open subset of the cotangent space to a 3-dimensional sphere (resp., 
a 3-dimcnsional two-sheeted hyperboloid) is symplectic. That diffcomorphism was 
discovered by Gyorgyi [9] (1968), re-discovered by Ligon and Schaaf [15] (1976) 
and discussed by Cushman and Duistcrmaat [5] (1997), who have shown that it is 
characterized by the three very natural properties: 

1. It maps the set of points in the phase space of the Kepler problem where the 
energy is negative (resp., positive) onto the tangent bundle of the 3-sphcrc 
(resp., the two-sheeted 3-dimensional revolution hyperboloid) with its zero 
section removed. 

2. It intertwines the Kepler and Delaunay vector fields (a rescaling of the geodesic 
vector field on the 3-sphere, or on the 3-dimensional hyperboloid). 

3. It intertwines the so(4)-momentum mappings of the Kepler and Delauney 
vector fields. 

We will see that the remarkable properties of that diffcomorphism appear as very 
natural consequences of the results presented in Section 2. 

We will also discuss the weak infinitesimal symmetries of the Kepler problem, 
and we will show that their set is a Lie algebroid with a zero anchor map, rather 
than a Lie algebra. 

After completion and submission of the present paper, the work of Gert Heckman 
and Tim de Laat [11], recently posted on arXiv, was indicated to us. The method 
used by these authors to explain the properties of the Gydrgyi-Ligon- Schaaf dif- 
fcomorphism rests on the same ideas as ours. We also learnt that conformally 
Hamiltonian vector fields were used, in the theory of bi-Hamiltonian vector fields, 
by A.J. Maciejewski, M. Prybylska and A.V. Tsiganov [16]. 

2. Conformally Hamiltonian vector fields. 

2.1. Notations and conventions. Unless another assumption is explicitly stated, 
all manifolds, functions, applications, vector fields and differential forms considered 
in this work will be assumed to be smooth, i. e differentiable of class C°° . 

Definition 2.1. Let A" be a vector field on a manifold M. The differential equation 
determined by X is the ordinary differential equation 

^ -*&*»• 

The flow of that equation is the map <&x , defined on an open subset Dx of R x M, 
with values in M, such that, for each x G M, the parametrized curve t >-> <&x{t,x) 
is the maximal solution <p of the above differential equation which takes the value 
x for t = 0. It means that 
d$(t x) 

— ±J—L = X($(t, x)) for each (t, x) G D x , *x(0, x) = x for each x G M , 
and that the open subset Dx of K x M is such that, for each x g M 

4 = {tgR| (t,x) g D x ) 
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is the largest open interval of R on which a solution ip of the differential equation 
determined by X satisfying </?(0) = x can be defined. 

Remark 1. The map <5>x is sometimes called partial flow to distinguish it from 
the full flow ^x, which is the map defined on an open subset of I x R x M, 
with values in M, such that, for each f el and x <G M, the parametrized curve 
t I— > ^x(t,to,x) is the maximal solution of the above differential equation which 
takes the value xq for t = to- The vector field X being time-independent, we have 
$x(t, 1 , xo) = &x(t — to,x ), so the full flow is determined by the partial flow. It 
is no longer true for ordinary differential equations determined by a time-dependent 
vector field; the full flow should then be used instead of the partial flow. 

2.2. Change of independent variable in ordinary differential equations. 

Lemma 2.2. On a manifold M, let X be a vector field and g a nowhere vanishing 
function. Consider the two ordinary differential equations 

d(p(t) 



dt 
#(s) 



g(<p(t))X(ip(t)), (*) 

x(m)- (**) 

Let us assume that there exists a smooth function a : R x M — > R such that for each 
solution ip of (*) 

-o-(t,<p(t))=g(<p(t)). (***) 
For each solution ip : I v — » M of (*) defined on the open interval I v o/R 7 let 

a v {t) = a(t,(f(t)) . 

Then the function is a diffeomorphism from I v onto another open interval o~,p(I v ) 
o/R ; and the map 

if) : -> M , s i-> i/j(s) = <p o a~ 1 {s) 

is a solution of the ordinary differential equation (**). 

Proof. The derivative of the function o~ v , at each poit t G I v , is gUpitf), which 
never vanishes. Therefore o~ v is a diffeomorphism. The chain rule shows that ip is 
a solution of (**). □ 

Remark 2. For a general smooth vector field X and a general smooth nowhere 
vanishing function g given on M, there may be no globally defined function a : 
I x M -} 1 verifying (***). In the Kepler problem, that function exists, as we will 
see in subsection 3.4, and is affine in the variable t. The following proposition gives 
some information about the existence of the map a. 

Proposition 1. On a manifold M , let X be a vector field and g : M — > R\{0} a 
nowhere vanishing function. We denote by &x '■ Dx M the flow of X (2.1). For 
each to £ K, there exists a smooth function a , defined on an open neighbourhood 
Wt of {t } x M mix M, such that for each solution ip : I v — x M of the differential 
equation determined by X defined on an open interval I v containing to, and each 
t G I v , 

j t o-(t,<p(t))=g(<p(t)). (*) 

The function a is not unique: any smooth function defined on M can be chosen for 
its restriction to {to} x AI . The function a is affine with respect to the variable 



4 



CHARLES-MICHEL MARLE 



t if and only if its restriction ot a : M — > R, o~t {x) = o~(to,x) is such that the 
Lie derivative with respect to the vector field X of the function g — a ta is constant 
along each integral curve of X . When the preceding condition is satisfied, the affine 
extension of a to the whole K x M, still denoted by a, satisfies (*) for all solutions 
if : I M of the differential equation determined by X and all t G /. 

Proof. Let us choose any to £ K. The map (9, x) M> (to — 6, x) is a diffcomorphism 
of R x M onto itself, which sends the open subset Dx on which the flow <&x is 
defined onto the subset W to = { (t, x) £ M X M \ (t - 1, x) £ D x }- Therefore W to 
is an open subset of R x M which contains {to} X M. Let at Q : M — > M be any 
smooth function. The formula 

a(t,x) = a to ° $x{to ~ t, x) + / jo$ x (t -t,x)dr 

Jto 

defines a smooth function a : Wt — > R whose restriction to {to} x M is (to, a;) i-> 
(to,at Q {x)) . For each (t,xo) € R x M such that (t — to,^o) € Dx, we have 

o-(t,<f>x(t-t (h x )) = cr to o $x(*o - - icb^o)) 



+ / 3o$ x (r-t,$ x (t-t ,x ))dT 

■/to 

o"t (a;o)+ / go$ x (T -t ,x )dT . 



Therefore 

— a(t,$(t-to,x )) = g{$x(t-t ,x )) . 

Since 1 i— > $x(t — to, 2^0) is the maximal solution of the differential equation deter- 
mined by X which takes the value Xq for t = to, we see that the map a satisfies 
condition (*) of the statement above. 

The map a is affine with respect to t if and only if its partial derivative with 
respect to t does not depend on t. We have 

da(t, x) , , _ . . ,, _ . N > 

g = -(do- to o$ x (t -t,x),Xo<f> x (t -t,x)) 



r d 

g{x) + J —(go<f> x ( T -t,x))dT. 



Taking into account 



we obtain 

da(t, x) 
dt 



d d 
— (go$x{T-t,x)) = - — (go<P x (T-t,x)) 



-(da to o $x(ta -t,x),Xo $ x (t -t,x)) 
+ g{x) - {{dg o $ x (r - t, x), Xo$ x ( T -t, x))\ 



T=t 

Setting 6 = to — t, we see that a is affine with respect to t if and only if, for all 

(9,x) e D x , 

(d(g - a ta )o$ x (0,x) 1 Xo $ x (0, a,)) = (d(g - a to )(x), X(x)) , 
which can also be written as 

C(X)(g - a to )($x(0,x)) = C(X)(g - a to )(x) , 
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where C(X) (g — at ) is the Lie derivative of the function g — ut a with respect to the 
vector field X. This equality expresses the fact that C(X)(g — (J t(l ) is constant along 
each integral curve of X . When the preceding condition is satisfied, the function 
a can be uniquely extended into a function, affine with respect to the variable t, 
defined on R x M, and one easily check that it satisfies condition (*) of the statement 
above for all solutions ip : I M of the differential equation determined by X and 
all tel. □ 

2.3. Hamiltonian and conformally Hamiltonian vector fields. 

Definition 2.3. Let (M,oj) be a smooth symplectic manifold and H : M — > R a 
smooth function. The unique vector field Xjj such that = — dH is called 

the Hamiltonian vector field associated to H, and H is called a Hamiltonian for Xh- 
Let g : M — > R\{0} be a nowhere vanishing function. The vector field Y = gXjj 
will be called a conformally Hamiltonian vector field, with H as Hamiltonian and 
<7 as conformal factor. The vector field V satisfies 

i(Y)u) = gi(X H )uj = -gdH . 

Theorem 2.4. Let (M, uj\) be a symplectic manifold, H : M — > R a smooth Hamil- 
tonian, X the associated Hamiltonian vector field. We assume that X is complete; 
in other words, its flow <&x is defined on the whole o/R x M. Let M° be an open 
dense subset of M , g : M R\{0} be a smooth, nowhere vanishing function and 
Y = gX be the conformally Hamiltonian vector field on M° , with Hamiltonian H 
and conformal factor g. Its flow will be denoted by $y. Let a : R x M° — > R be a 
smooth function such that for each maximal solution Lp of the differential equation 
determined by Y , 

— it — =3im- 

We assume that there exists on M° another symplectic form u>2 such that 

i{Y)uj 2 = -dH . 

In other words, the vector field Y is both Hamiltonian with respect to loi with H as 
Hamiltonian and conformally Hamiltonian with respect to u>\ with the same H as 
Hamiltonian and with g as conformal factor. 
Under these assumptions, the map 

H:M°^M, x i-> E(x) = $x(-ct(0,x),x) 

is a symplectic diffeomorphism from (M ,W2) onto an open subset of (M,uii), equi- 
variant with respect to the flow of Y on M° and the flow of X on M , that is 

and for each (t,x) in the open subset o/R x M° on which $y is defined 

$ x (t,E(x))=E(<f> Y (t,x)). 

Proof. The maximal integral curve of the differential equation determined by Y, 
which takes the value xq for t = to, is 

t i y $y(t - t ,Xo) . 

The same geometric curve in M , parametrized by s = a(t,x) instead of t, is an 
integral curve of the differential equation determined by the vector field X. The 
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values of the parameter s which correspond to (to, xo) and to (t, $(i — to, Xo)) are, 
respectively, 

so = <j(t ,x ) and s = a(t,$(t - t ,x )) . 
Since $x is the flow of the vector field X, we have 

$ y (t - t , x ) = $x (v(t, - t , x a )) - a(t , x ) , £0) ■ 

Let T:lx M° — > K x M be the map 

(t,x) ^ T(t,x) = (t,$ x (t-a(t,x),x)^) . 

We are going to prove that for all t, to, xo such that 5>y (t — t ,xo) is defined, 

r(t,® Y (t-t ,xo)) = (t,$x{t-<r(to,xo),xo)) ■ 

The above formula expresses the fact that T maps the graph of the integral curve 
t i — ^ — io,£o) of the vector field Y which takes the value xo for t = to, into 

the graph of the integral curve t H> $>x(t — a(to,xo),xo) which takes the value Xq 
for t = <r(t ,xo)- Observe that T associates, to the point of the integral curve of 
Y reached for the value t of the parameter, the point of the integral curve of X 
reached for the same value s = t of the parameter. 

Replacing x by & Y (t — to, xq) in the formula which defines T, we get 

r(t,$Y(t-t ,X )) = ^t,$x(t-<7(t,<I> Y {t-to,Xo)),<l>Y{t-to,Xo)^ . 

But we have shown that 

$ Y {t-to, x ) = $x (v(t, $ Y (t - t ,x )) -a(t ,x ), x ^j ■ 

Therefore 



(t - a(t, <P Y (t - t , x )),$ Y (t - t Q , x )) 

= ®x (t-v{t,®Y{t-to,x )),$ x (a(t,$ Y (t-t ,Xo)) -a(t ,xo), xo)^ 

= <!>x(t- u(to,xo),xo) ■ 
We have proven that T(i, $(£ — to,xo)) = [t, ty(t — er(to, xo), xo) ). Since Y is a 



Hamiltonian vector field on (M°,W2), with H as Hamiltonian, the kernel of the 
closed 2-form on R x M° 

Ld2 = (-^2 — dH A dt 

is the rank-one completely integrable distribution generated by the nowhere van- 
ishing vector field 

d 

ot +Y > 

where t denotes the coordinate function on the factor R. We have indeed 
^ + Yj (uj 2 - dH Adt) = dH - dH = . 

Similarly, since X is a Hamiltonian vector field on (M, uji), with H as Hamiltonian, 
the kernel of the closed 2-form on R x M 

uji = uji — dH A ds 
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is the rank-one completely integrable distribution generated by the nowhere van- 
ishing vector field 

d 

7T + X > 
os 

where the coordinate function on the factor R is now denoted by s. We have indeed 
d \ 

— + X - dH A ds) = dH - dH = . 
os J 

We recall that T:lx AI° — >• R x M maps injectively each leaf of the foliation of 
K x M° into a leaf of the foliation off xM. On the manifold R x M°, the 2-forms 
UJ2 and T*wi both have the same kernel, since their kernels determine the same 
foliation. Each of these 2-forms is therefore the product of the other by a nowhere 
vanishing function. This function is in fact the constant 1, because 

T*s = t, ,T*H = H, so T* (dH A ds) = dH A dt , 

and we have Y*wi = cJJ. 

Restricted to {0} x M°, the map T becomes 

(0,o:) i — ^ T(0,o:) = (o,^(- a(0,x),x)) = (0,S(x)) . 

Since u>2 and lo\ are the forms induced, respectively, by uji on {0} x M° and by ui\ 
on {0} x M, we have 

S*Wi = UJ2 ■ 

The 2-forms uj\ and L02 being nondegenerate, that proves that the map S is open. 
But its geometric interpretation proves that 3 is injective: for each xo € A/ , 
H(a;o) is the point of M reached, for the value s — of the parameter s, by the 
integral curve ip of X whose value for s = s = cr(0,a;o) is i/j(sq) = x ; it always 
exists, because X is assumed to be complete. Being open and injective, S is a 
diffcomorphism of M° onto an open subset of M, and more precisely, since S*wi = 
0J2, a symplectic diffcomorphism of (M , uj%) onto an open subset of (M,lui). 

We have seen that T is equivariant with respect to the flows of the vector fields 
d d 

— + Y on R X M° and — — h X on R x M. By projection on the second factor, we 
ot os 

see that 2 is equivariant with respect to the flows of Y on M° and X on M. □ 
3. Application to the Kepler problem. 

3.1. The equations of motion of the Kepler problem. In the physical space, 
mathematically modelled by an Euclidean affine 3-dimensional space £ , with asso- 
ciated Euclidean vector space c , let P be a material point of mass m subjected 
to the gravitational field created by an attractive centre O. The physical time is 
mathematically modelled by a real independent variable t. Let 



d~^ 

r = 1 1 r' 1 1 ; p = rr, 

The force exerted on P i 



~F" = OP" ; r = |~F"j ; = m : p = II 

dt 



7-- 



km 



where fc is the constant which characterizes the acceleration field created by O. 
The equations of motion are 

cTF* jT* dp km~r^ 
dt m dt r 3 
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The second equation above becomes singular when r = 0, so we will assume that 
the massive point P evolves in £\{0}. That space will be called the configuration 
space of the Kepler problem. 

3.2. First integrals. The above equations of motion are those determined by a 
Hamiltonian, time-independent vector field on a symplectic 6-dimensional manifold, 
called the phase space of the Kepler problem: it is the cotangent bundle T* (£\{0}) 
to the configuration space. We will identify the tangent and cotangent bundles by 
means of the Euclidean scalar product. The Hamiltonian, whose physical meaning 
is the total energy (kinetic plus potential) of P, is 

„ . • s > . P 2 Tnk 

E = E(?, = £ • 

Am r 

The energy E is a first integral of the motion, i.e., it is constant along each integral 
curve. 

The group SO(3) acts on the configuration space £\{0} by rotations around O; 
the canonical lift of that action to the cotangent bundle leaves invariant the Liouville 
1-form and its exterior derivative, the canonical symplectic form of T* (£\{0}) . It 
also leaves invariant the Hamiltonian E. The lifted SO(3)-action is Hamiltonian, 
and the corresponding momentum map is an so(3) "-valued first integral. Since the 
Euclidean vector space £ is of dimension 3, once an orientation of that space has 
been chosen, we can identify the Lie algebra so(3) and its dual space so(3)* with the 
vector space ~t itself. With that identification the bracket in the Lie algebra so(3) 
becomes the vector product, denoted by (~u , v) i— > u x v\ the coupling between 
spaces in duality becomes the scalar product (ll , v) i— > u ■ v ; the momentum map 
of the Hamiltonian action of SO (3) is (up to a sign change) the well known angular 
momentum 

L = ~F* x JI . 

In addition to the energy E and the angular momentum L , the equations of motion 
of the Kepler problem have as a first integral the eccentricity vector e discovered 
by Jakob Herman (1678-1753) [12, 2], improperly called the Laplace vector, or the 
Runge-Lenz vector, whose origin remained mysterious for a long time: 

*__!* + t42_fJt-iW-2£*. 



r m 2 k \m 2 k r J m 2 k 

Using the well-known fact (the first Kepler's law) that the orbit (the curve described 
by P as a function of time) is a conic section with O as one of its foci, the eccentricity 
vector has a very simple geometric meaning: it is a dimensionless vector, parallel 
to the major axis of the orbit, directed from the attracting centre towards the 
perihelion, of length numerically equal to the eccentricity of the orbit. 

3.3. The hodograph. We will not recall in full how the equations of motion of the 
Kepler problem can be solved, since it is done in several excellent texts [1, 4, 8, 18]. 
Let us however briefly indicate the proof, due to William Rowan Hamilton (1805- 
1865) [10], of an important fact: the hodograph 1 of each motion of the Kepler 
problem is a circle or an arc of a circle. 



The hodograph of the motion of a particle moving in an affine space E is the curve, drawn in 
the associated vector space £ , by the velocity vector it of the particle, as a function of time. 
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Let us first look at solutions for which at the initial time to, and ~jt are not 
collinear, so L (to) ^ 0. Since L is time-independent, 1^ and ~p* are never collinear. 
We choose an orthonormal frame, positively oriented, with O as origin and unit 
vectors e^, and et with L parallel to ej. We have 

L = Let , where L is a constant. 

The vectors et p remain for all times parallel to the plane xOy. Let 9 be the 
polar angle made by with e x . We have 

r = r cos + r sin , 

-> /dr „ d0 \ fdr . n d9 \ 

p = m — cos v — r — sin 9 ) el + m \ — sin # + r — cos w e„ 

-> 2 d(9^ 
1/ = mr — e, . 
dt 

Therefore 

2 d6' „ 
L = mr — = Constant . 
dt 

This is the second Kepler law, also called law of areas, since is the area swept 

. . . 2m 

by the straight segment OP during an unit time, with sign + if 9 increases with 
time and — if it decreases. Observe that 9's variation is strictly monotonic. So we 
can take 9 as independent variable, instead of time t. We may write 



dp d~f dt m 2 k . , . 

-dJ = lttd9 = - — ( cos ^ + sm ^) 



This ordinary differential equation for the unknown ~p* , which no more involves r , 
can be readily integrated: 

— ± m 2 k . . — ^ — ). — y 
p = ——(— sm0e x + cos9e y ) + c , 
Li 

where c is a (vector) integrating constant. We will choose such that 

C C€y , 

where c is a numeric constant of the same sign as L. 

With O as origin let us draw two vectors in the plane xOy, the first one (constant) 
being equal to ~c , and the second one (which varies with 9) equal to ~p*. The end 
point of that second vector moves on a circle whose centre is the end point of the 
vector equal to c , and whose radius is 

m 2 k 

K = W 

This circle (or part of a circle) is (up to multiplication by m) the hodograph of the 
Kepler problem. 

A short calculation leads to the following very simple relation between the energy 
E of a motion, the radius 71 of its hodograph and the distance |c| from the attracting 
centre O to the centre of the hodograph: 

2mE = c 2 - K 2 . 
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Observe that the right-hand side c 2 — 1Z 2 is the power 2 of O with respect to the 
hodograph. 

3.4. The Levi-Civita parameter. Let a be the function, defined on the product 
with R of the phase space of the Kepler problem, 

mk 

A short calculation using the equations of motion shows that for any solution t i— > 
(r(t],p(t\) of the Kepler problem, 

da-(r^j,p(tj) _ 1 

dt ~~ r(i) ' 

With s(t) = a(r(t),p(t)) as the new independent variable, instead of the time t, the 
equations of motion become 




ds m ds r 2 {s) 



I will call s the Levi-Civita parameter: it was introduced, as far as I know for the 
first time, by Tullio Levi-Civita in [13], by the differential relation which expresses 
ds as a fuction of dt and r{t). The integrated formula which gives o~(t, r, p) is in 
the paper [21] by Jean-Marie Souriau, and in Exercise 8, chapter II of the book by 
R. Cushman and L. Bates [4]. This formula has probably been known before for a 
long time in the Celestial Mechanics community, but I do not know who found it 
for the first time. 

With the Levi-Civita parameter as new independent variable, the system is no 
longer Hamiltonian, but rather conformally Hamiltonian. The Levi-Civita param- 
eter will be used in subsection 3.14 in a slightly different context: we will define a 
diffcomorphism S from the phase space of the Kepler problem restricted to negative 
(resp. positive) values of the energy, onto an open dense subset of the cotangent 
bundle to a three-dimensional sphere (resp., to one sheet of a two-sheeted three- 
dimensional hypcrboloid) . On this new phase space equipped with its canonical 
symplectic form, the image of the vector field which determines the equations of 
motion of the Kepler problem will be conformally Hamiltonian instead of Hamil- 
tonian, while the image of the vector field transformed by the use of the Levi-Civita 
parameter as independent variable will be Hamiltonian. 

3.5. The Gyorgyi-Ligon-Schaaf symplectic diffeomorphism. G. Gyorgyi [9] 
gave the expression of a symplectic diffeomorphism from the phase space of the 
Kepler problem restricted to the negative (resp. positive) values of the energy, onto 
an open subset of the cotangent bundle to a three-dimensional sphere (resp. two- 
sheeted hyperboloid). He obtained this diffeomorphism in two steps. First, following 
the ideas of Fock [7], employed two years later by Moser in his well known paper [19], 
he arranged the components of the first integrals E, L and of the Kepler problem 
into a 4 x 4 matrix, and used that matrix to define two 4-dimensional vectors ~f} 
and 7^ (formulae 2.15 and 2.16 of Gyorgyi's 1968 paper [9]). I believe that the map 

2 In plane Euclidean geometry, the power of a point O with respect to a circle C is the real 
number oX.O^i, where A and B are the two intersection points of C with a straight line T> through 
O. That number docs not depend on T> and is equal to ||0(3|| 2 — 1Z 2 , where C is the centre and 
7Z the radius of C. 
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(r,p) i-> (p,ft) obtained by Gyorgyi after difficult to follow calculations, is the 
map, maybe rescaled, I have called S 1-1 in Subsection 3.f 2 below. Second, in Section 
entiteled "Getting back the time t; Bacry's generators" of his 1968 paper, Gyorgyi 
composes the map ( r , p) i— > (p, vr ) with a suitably chosen map deduced from the 
flow of the Kepler vector field. He did not prove that the map he finally obtained 
is symplectic. This result was proven by Ligon and Schaaf [15] who, eight years 
later, rediscovered the same diffcomorphism. Twenty years later, R. H. Cushman 
and J. J. Duistermaat [5] discussed the properties of that diffeomorphism, gave new 
proofs of its remarkable properties and stated the nice characterization indicated in 
the Introduction. 

For negative values of the energy, the Gyorgyi-Ligon-Schaaf diffcomorphism is 
the map which, to each element (r,p) of the phase space of the Kepler problem 
such that E(r ,p) < 0, associates the element ((£o, £ ), {Wo,y)) of R 4 x R 4 



y/-2mE(?,J) ^ . frp 2 \ 

£o = ; r.p sin ip + — - - 1 cos ip , 

mk \mk J 

-> . , j-2mE{T,J) _+ 

£ = I r~ P sin w H ; r p cos w , 

\ r mk J mk 

—> —> ( rp 2 \ . 

r/ = - r . p cos tp + —======= — - 1 \ smtp , 

^-2mE(T,p) \mk 



mk 

^-2mE{y, f) \ r mk 



if = ^^=^^^=^== ( — \ cos (p + r~0 sin ip . 



In these formulae tp is the angle given, as a function of (~r^, ~$), by 

_> y/-2mE(T, J) 

p > = r=r^ ■ 

mk r . p 

The quantities ~£ = (£1,62, £3) an(1 ^ = (^1 , ? ?2, arc vectors of R 3 , while 
(Co, £ ) = (Co, £1,6, &) and (i]o,~rf) = (rjo, f?i, ^2,^3) are vectors of R 4 , such that 

£o 2 + lltll 2 =£o 2 + £? + £ 2 2 + £3 2 = i ) 
-> — > 

£0% + £•?/ = &'7o + £im + 6?72 + £3% = , 

V 2 o + \\^f\\ 2 = Vo+vi + V 2 2 +vi >0. 

In other words, (£o, £ ) is a point of the sphere S* 3 of radius 1 centered at the origin, 
embedded in R 4 (endowed with its usual Euclidean metric), and (??o,^f) is a non- 
zero vector tangent to that sphere at point (£o, £ ). Using the usual scalar product 
of R 4 to identify vectors and covectors, we may consider ((£ , £ ), (rjOiV)) as an 
element of the cotangent bundle to the sphere S 3 minus its zero section. 

For positive values of the energy, the Gyorgyi-Ligon-Schaaf diffcomorphism is 
given by similar formulae, the trigonometric functions sin and cos being replaced 
by the hyperbolic functions sinh and cosh and the Euclidean metric on R 4 by the 
Lorentz pseudo-Euclidean metric. 

We will show in what follows that the Gyorgyi-Ligon-Schaaf symplectic diffco- 
morphism can be easily obtained, and that all its properties can be proven, by an 
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adaptation of Moscr's method for the regularization of the Kepler problem and ap- 
plication of Theorem 2.4. That method rests on the fact that hodographs of the 
Kepler problem are circles or arcs of circles, and that the stercographic projection 
maps circles into circles. We do not need to know in advance that the eccentricity 
vector is a first integral of the Kepler problem, since this property is an easy 
consequence of Moser's method. 

3.6. Stereographic projection. In 1935, V. A. Fock [7] used an inverse stere- 
ographic projection for the determination of the energy levels of hydrogen atoms 
in quantum mechanics. More recently J. Moser [19] used the same idea for the 
regularization of the Kepler problem. His method yields a correspondence between 
oriented hodographs of motions of the Kepler problem with a fixed, negative en- 
ergy, and great oriented circles of a three-dimensional sphere. As explained by 
D. V. Anosov [1] and J. Milnor [18], the method can easily be adapted for motions 
with a fixed, positive energy: instead of a three-dimensional sphere, one should use 
one sheet of a two-sheeted three-dimensional hypcrboloid of revolution, and one 
obtains a correspondence between oriented hodographs of motions of the Kepler 
problem with the chosen positive energy, and connected components of great hy- 
perbolas drawn on that sheet of hyperboloid (that means intersections of that sheet 
of hypcrboloid with a plane containing the symmetry centre) . 

To deal with cases E < and E > simultaneously, we introduce the auxiliary 
quantity 




1 if £<0, 
-1 if E>0. 



Let (O, eZ, iy , el) be an orthonormal frame of £. with O as origin. We add to the 

basis (et,et),et) of ~~£* a fourth vector e~l and we denote by h the corresponding 
coordinate. So we get a 4-dimensional affine space J- ' . The physical space £ will be 
identified with the affine subspace of JF determined by the equation h = 0. 
In J 7 , for each p > 0, let Q p be the quadric defined by 

h 2 + C(.* 2 + y 2 + z 2 ) = p 2 . 

It is a sphere if £ — 1, a two-sheeted hyperboloid if C = — 1- We will see that there 
is a value of p particularly suited for each value of the energy E. 

Let N be the point with coordinates (x = y = z = 0, h = p). The stereographic 
projection (usual if ( = 1, generalized if £ = —1) of the quadric Q p minus point 
TV on the space £ is the map which associates, with each point M € Q P \{N}, 
the intersection point m of the straight line which joins N and M, with £. See 
Figure 1. If £ = 1 that map is a diffeomorphism from Q p \{N} onto £. If £ = — 1 
it is a diffeomorphism from Q p \{N} onto the open subset of £ complementary to 
the 2-sphere of centre O and radius p. The upper sheet (h > 0) of the hyperboloid 
(minus point A'') is mapped onto the outside of that sphere and the lower sheet 
(h < 0) onto its inside. 

3.7. Cotangent lift. Since the (maybe generalized) stereographic projection is a 
diffeomorphism, it can be uniquely lifted to the cotangent bundles in such a way 
that the pull-back of the Liouville form on T*£ is equal to the Liouville form on 
T*(Q p \{N}). We apply this construction with Om = p ', the linear momentum of 
a Keplerian motion of energy E ^ 0. 
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Figure 1. Stcrcographic projection and its generalization 



Each 1-form on £ can be written • dOm, where is a vector field on £, or, 
since Om = p, 

f ■ d~f = r x dp x + r y dp y + r z dp z . 
Recall that ( = 1 if E < and ( = — 1 if E > 0. The equation of Q p is 

otl ■ OaI = x 2 + y 2 + z 2 + Qh 2 = Cp 2 . 
Each 1-form on Q p can be written 

■ dOM = W x dx + W v dy + W z dz + QW h dh , 

where is a vector field defined on Q pi tangent to that quadric. Therefore, 

■ OM = W x x + W v y + W z z + (W h h = . 

The cotangent lift of the (maybe generalized) stereographic projection will be 
denoted by S p . It maps each pair (M, W) made by a point M £ Q P \{N} and a 
vector W tangent to Q p at that point, onto a pair of two vectors of £ , in 

such a way that 

■ dotl = ~? ■ d~f . 

The formulae which give (~]} , as functions of (OM, are easily obtained after 
some calculations: 

~f = —^—oti, with h 2 + aojif = p 2 , 
p-h 



where we have set 



OM = Op + hel, W = W 3 + W h el 



— i 
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The formulae for the inverse transform are: 



Op 



p 2 + (p 2 

p 2 ~Cp 2 



h = p p 2 + Cp 2 : with p= = V p * +p y +p2 ^ 



P" + c,p~ c, r ■ p -4 
w 3 = -^-?- — P 



These formulae prove that if pt = Omi and p| = Om-i are two collinear vectors in 

E such that 

Om\.Om,i = — (p 2 , 

their images, by the inverse stercographic projection, are two points Mi and M2 
symmetric of each other with respect to O. That property may also be proven with- 
out calculation as a consequence of properties of the transformation by inversion. 

Wc recall that if a straight line through O meets the hodograph of a Kcplerian 
motion at two points m 1 and m 2 , Om\.Omi is the power (see the second footnote 
in Subsection 3.3) of O with respect to the hodograph, and is equal to 2mE, where 
E is the energy of the motion. Therefore the above results shows that the inverse 

(p 2 

streographic projection maps the hodographs of motions of energy on curves, 

2m 

drawn on the quadric Q p , which are symmetric with respect to O. These curves are 
great circles of the sphere Q p if ( = 1, and great hyperbolas of the hyperboloid Q p 
if £ = — 1. The name "great hyperbola" is used by analogy with "great circle": it 
means the intersections of Q p with a plane containing the symmetry centre O. As 
a consequence, the cotangent lift S p x of the inverse stereographic projection maps 

Cp 2 

the hodographs of motions of the Kepler problem of energy E = on geodesies 

2m 

of the quadric Q p . The hodographs of motions with a different value of the energy 
are mapped on circles (if £ = 1) or hyperbolas (if C = — 1) drawn on the quadric 
Qp, but which are not great circles, or great hyperbolas. 

So we see that for each value of the energy E, there is a particularly well suited 

Cp 2 

value of p: it is the unique p > such that E = , with C = 1 if E < and 

2m 

C = -1 if E > 0. 

Remark 3. We have defined the mapping S p : T*(Q p \{N}) -> T*£ so that 

(S p )*C? -d-f) = $-dOl3. 

Since the symplcctic form of T*E is d(~f ■ drf) and that of T*(Q P \{N}) d{\& ■ 

dOM), the map S p is an anti-symplectic diffeomorphism, rather than a symplectic 
diffcomorphism. 

The transformed Hamiltonian (S P )*E = E o S p is 

EoS p = ^ + C -^(\\^\\- k 4 

2m m (p-h)\\W\\ \ P 

where we have set 

ll^ll = sjw% + W 2 + W 2 + (W 2 . 
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The expression of Eo S„ shows that E (T^, = if and only if S„ 1 (~r^, = 

2m 1 

/ y y kui^ 

[OM, Wj is such that \\W\\ = — r . By differentiation we obtain 



P- 

d (E o Sp) = , ^ d (ll^ll) + ( - ^ 



m(p-h)\\W\\ y 7 V P / \m(p-/i)||W , ||/ 

CP 2 

For a motion of energy E = , the second term of the right hand side of the 

2m 

/ l kTTl \ 

above equality vanishes, since ( IIW1I tt- I = 0, and we see that S^ 1 transforms 

V p 2 J 

the Hamiltonian vector field of the Kepler problem into a conformally Hamiltonian 
vector field on T*(Q p \{N}) minus the zero section, with Hamiltonian ||W|| and 

-Cp 3 

conformal factor We introduced a minus sign in the conformal 

m{p-h)\\W\\ 

factor to account for the fact that S p 1 is anti-symplectic (see Remark 3). 

3.8. Motions of zero energy. For Keplerian motions of energy zero the stere- 
ographic projection should be replaced by an inversion with pole O and ratio I. 
Formulae for its cotangent lift arc 

t= =L oti, 

l|OM|| 2 

\\otl\\ 2 -> 2V^ -otl \ 

7* = " " W OM . 

The correspondence (jf, M> (OM, W) being involutive, the formulae for the 
inverse transformation are 

Denoting by Sq the cotangent lift of the inversion of pole O and ratio I, we easily 
obtain the transformed Hamiltonian 

p c t 2m2k 

E O 5n — =5 r \\W\\ - 



2m\\W\\ \\OMf 
Now W is a vector of £ , with 3 components, so we have set 



\\w\\ = ^w* + w* + w*. 

As above this result proves that E(r,p) = if and only if Sq 1 (r, p) = 

(Olft, is such that \\V\f\\ = 
By differentiation we obtain 

d 



d(EoS ) = _/ = d (0\\) + 

2m\\W\\ \\OM\\ 2 V 1 \ 1 ) \2m\\W\\ \\OM\\ 2 , 
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For a motion of energy E = 0, the second term of the right hand side of the above 

i 2fcm 2 \ 

|| Ty || — J = 0, and we see that Sq 1 transforms the 

Hamiltonian vector field of the Kepler problem into a conformally Hamiltonian vec- 
tor field on T* (£\{0}) minus the zero section, with Hamiltonian \\W\\ and conformal 
-I 2 

factor =r — r — . We introduced a minus sign in the conformal factor to ac- 

2m||VK|| ||OM|| 2 

count for the fact that Sq is anti-symplcctic (Remark 3 applies to So as well as to 

s P ). 

3.9. Infinitesimal symmetries. Let us first recall some properties about infini- 
tesimal symmetries. Consider a vector field X on a manifold M. Another vector 
field Z on M is said to be an infinitesimal symmetry of the differential equation de- 
termined by X (or, in short, of X) if [Z, X] = 0; it is said to be a weak infinitesimal 
symmetry of X if at each point x <E M, [Z,X](x) and X(x) are collinear. If Z is 
an infinitesimal symmetry or a weak infinitesimal symmetry of X and h a smooth 
function, Z is a weak infinitesimal symmetry of hX , since 

[Z, hX] = h[Z, X] + (C(Z)h) X , 

whichs shows that at each point x £ M, [Z,hX](x) and X(x) are collinear. When 
X = Xh is a Hamiltonian vector field on a symplcctic manifold (M, uj), with a func- 
tion H as Hamiltonian, a vector field Z on M such that C{Z)us = and C(Z)H = 
is an infinitesimal symmetry of Xh and a weak infinitesimal symmetry of hXjj 
for any smooth function h. When in addition Z is Hamiltonian, its Hamiltonian 
is a first integral of Xh and of gXn- Conversely, if / and g are two smooth first 
integrals of the Hamiltonian vector field Xh, the vector field gXf, Xf being the 
Hamiltonian vector field with Hamiltonian /, is an infinitesimal symmetry of Xh 
and a weak infinitesimal symmetry of hXn for any smooth function h. We have 
indeed 

[ 9 X f ,hX H } = gh[X f , X H ] + g{i(X f )h) X H - h(i(X H )g) X f 

= ghX {f , H} + g{f, h} X H - h{H, g} X f = g{f, h} X H , 

since {f,H} = {H,g} = 0. 

For the vector field which determines the equations of motion of the Kepler prob- 
lem, we already know a three-dimensional vector space of infinitesimal symmetries: 
it is made by the canonical lifts to the cotangent bundle of infinitesimal rotations 
of the configuration space £\{0} around the attractive centre O. These canonical 
lifts are Hamiltonian vector fields, and their Hamiltonians (linear combinations of 
the components of the angular momentum L) are first integrals of the equations 
of motion. The results obtained in Subsections 3.7 and 3.8 show that other weak 
infinitesimal symmetries exist. In these subsections, for each value e of the energy 
E, we have built a symplectic diffeomorphism 3 of the phase space of the Kepler 
problem onto an open subset of another symplectic manifold: the cotangent space 
to a three-dimensional sphere when e < 0, the cotangent space to £ when e = and 
the cotangent space to a three-dimensional hyperboloid when e > 0. That sphere, 
affine space £ or hyperboloid will be called the new configuration space, and its 
cotangent bundle (with the zero section removed) the new phase space. 



3 Or rather an anti-symplectic diffeomerphism, but by a change of sign it can be transformed 
into a symplectic diffeomorphism 
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For each energy level e, the new phase space has a 6-dimensional Lie group of 
global symmetries: SO(4) when e < 0, the group SE(3) of orientation-preserving 
isomctries (rotations and translations) of £ when e = and the Lorentz group 
SO(3, 1) when e > 0. The direct image of the Hamiltonian vector field of the Ke- 
pler problem, restricted to each energy level of the phase space, is equal to the 
restriction, to the corresponding submanifold of the new phase space, of a confor- 
mally Hamiltonian vector field, whose Hamiltonian is invariant by the action of 
the 6-dimensional Lie group of symmetries of the new phase space. Therefore, on 
each energy level of the phase space of the Kepler problem, six linearly independent 
weak infinitesimal symmetries exist, of which three are the already known linearly 
independent infinitesimal symmetries associated with the cotangent lifts of infini- 
tesimal rotations of £ around three non-coplanar axes through O. The other three 
correspond to the additional symmetries of the considered energy level of the new 
phase space. We will see that they depend smoothly on the value of the energy. 
Therefore, on the phase space of the Kepler problem considered as a whole, there 
exist 6 linearly independent vector fields, three of them being infinitesimal symme- 
tries and the other three weak infinitesimal symmetries. We will see that these six 
vector fields, which of course are tangent to each energy level, arc Hamiltonian and 
are all infinitesimal symmetries of the Kepler vector field. So their Hamiltonians 
are first integrals of the equations of motion. 

In Subsection 3.10 we will see that these infinitesimal symmetries do not span 
a Lie algebra, but rather a fibered space in Lie algebras [6] over the base R, the 
value of the energy running over that base, in other words a Lie algebroid [3, 17] of 
a special kind, with a zero anchor map. 

Let us recall [14, 20] that if $ : G x M — > M is a left action of a Lie group G 
on a manifold M, its canonical lift to the cotangent bundle is a Hamiltonian action 
$ : G x T*M -> T*M with an equivariant momentum map J : T*M — s- Q* given 



In this formula, 7Tm : T*M — > M is the canonical projection; X G Q, the Lie algebra 
of G; £ G Q* , the dual space of Q\ and X M is the fundamental vector field on M 
associated to X, defined by 



The momentum map J is a first integral of any Hamiltonian vector field Xh or con- 
formally Hamiltonian vector field gXji on T*M whose Hamiltonian H is invariant 
under For each X G Q the Hamiltonian vector field on T*M whose Hamiltonian 
is £ i — y (J(£),X) is an infinitesimal symmetry of Xjj- 

We apply these results, with for 4> the action on the new phase space of its 
symmetry group. As a vector space, the Lie algebra of the symmetry group is 
canonically isomorphic to but with a Lie algebra bracket [ , ] e which 

depends on the energy level e. That bracket will be determined in Subsection 3.10. 
Two kinds of infinitesimal symmetries of the new phase space exist. The first kind 
is made by the fundamental vector fields which correspond to elements (ui,0) of 
These vector fields are the canonical lifts to the cotangent bundle of the 
infinitesimal rotation of £ around the axis through O parallel to u{. The second 



by 



(J(0,X) = (t,X M (ir M (0))- 




for each x G M . 
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kind is made by the fundamental vector fields which correspond to elements (0, u?) 
of ~? x~t. These vector fields arc the canonical lifts to the cotangent bundle of the 
following infinitesimal transformations of the new configuration space: 

• when e < 0, the restriction to the sphere Q p , with p 2 = — 2me, of the infini- 
tesimal rotation of x (Ixel) in which the plane spanned by U2 and 

rotates, while the vector subspace of ~~£* orthogonal to «2 remains fixed; 

• when e = 0, the infinitesimal translation of £ parallel to U2', 

• when e > 0, the restriction to the hyperboloid Q p , with p 2 = 2me, of the 
infinitesimal Lorcntz transformation of x (K x eh) in which the plane 
spanned by U2 and et is transformed, while the vector subspace of ~£* orthog- 
onal to Tia remains fixed. 

The space ?x^ will be identified with its dual space by using the scalar product 
on each factor. The momentum map, which will be denoted by K e , can be written 
K e — (Kie, K2e), each component taking its values in £ , 

In the three cases e < 0, e = and e > 0, the Hamiltonian H has the same 
expression 

(oaI,#) ^ #(oaI,w^) = \\vfr\\ . 

It is equivalent to the classical Hamiltonian — - — of a particle moving freely on 

Q p when e ^ or on £ when e = 0. It is indeed invariant under the action $. 
The first component K\ e of the momentum map is 



|0/j x W3 when e^0, 



[OMxW when e = 0, 
and its second component is 

U{hW^-W h Op) when e ^ 0, 
1 W when e = 0. 



As above, we have set when e ^ 



OM = Op + het, = W 3 + W h e^ 



ii 



Let us now come back to the phase space of the Kepler problem. The composed map 
T — K~ n S*^ 1 , which has two components J\ e = K\ e o S" 1 and J 2e = K 2e Sp 1 



is the momentum map o f the Lie algebra action of infinitesimal symmetries. Using 
the formulae which give OA? and W as functions of 1^ and jt, we obtain for the 
first component, in the three cases e < 0, e = and e > 

Ji e = ~0 x "7* . 

Its expression does not depend on e. Up to a change of sign, it is the angular 
momentum L . For the second component we obtain 

when e 7^ 0, 
when e = 0. 
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Using the value of the energy E = — and, when e / 0, the equality 

2m r 

(p 2 = —2mE, these formulae become 

: £ when e^O, 



J 2 e = I fW* 

— — £ when e = 0, 

where is the eccentricity vector. Since the energy EC? , ~p ) is a first integral, we 
see that by choosing I = 2m 2 k and, for e 7^ 0, by multiplying the weak infinitesimal 

, . 1 ,1 , ,. y/-2tmE(T,T) 

symmetries by the smooth function — , we can arrange things so 

m^k 

that the second component of the momentum map becomes 

J2e = £ j 

which no longer depends on the energy level e. Therefore, although we have defined 
them separately on each energy level, the three additional weak infinitesimal sym- 
metries smoothly depend on the energy level, since they are the Hamiltonian vector 
fields whose Hamiltonians are the components of £ ; they are true (not only weak) 
infinitesimas symmetries since they are Hamiltonian vector fields whose Hamiltoni- 
ans are first integrals of the Kepler equations of motion. 

3.f0. The energy-momentum space and map. The Poisson brackets of the 
components of the angular momentum L and of the eccentricity vector it in an 



orthonormal frame (e£, e^,et) of ~~t 



are 



{L X , Ly} 


= -£., 


{Ly, L z } 


— — L x , 


{L Z ,L X } = 


L y ; 


{L x ,£ x } 


= 0, 


{L X , £y} 


= Sz , 


{L x ,e z } = 


£y ; 


{L y , e x } 


= £z , 


{Ly,£y} 


= 0, 


{L y ,z z } = 




{L z ,£ x } 




{ L z,£y} 




{L z ,e z } = 


0; 




- 2E L 
m 3 k 2 ' 


{ e yj e z} 


- 2E L 
m 3 k 2 ' 


{£z,e x } = 


2E L 

m 3 k 2 y 



The family of fuctions on the phase space of the Kepler problem spanned by the 
components of L and of £ is not a Lie algebra, because the function E appears 
in the right hand sides of equalities in the last line. Of course, if we replace £ 

by ; we get a Lie algebra isomorphic to so(4) when E < and to so(3, 1 ) 

when E > 0. But in our opinion, this is not a good idea, because if we do this the 
energy level E = is lost and the true geometric nature of the family of infinitesimal 
symmetries, on the whole phase space of the Kepler problem, is hidden. The above 
formulae make up the bracket table of smooth functions defined on the dual of a 
Lie algcbroid (with zero anchor map), whose base is R (spanned by the coordinate 
E) and standard fibre R 6 = R 3 x M 3 (spanned by the components of L and e as 
coordinate functions). 

The energy-momentum map is the map J , defined on the phase space of the 
Kepler problem, with values in R x £ x £ , 

H- JC?,^) = (Erf, (?,?),-?(?,?)) . 
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The space R x will be called the energy-momentum space. As seen above, 

that space is the dual of a Lie algebroid; it should be considered as fibered over its 
first factor R, the fibre {e} x txt over each point e G R being equipped with 
a linear Poisson structure which smoothly depends on e. With this structure, that 
fibre is the dual space of a Lie algebra which, as a vector space, can be identified 
with fx?, with a Lie algebra bracket [ , ] e which smoothly depends on e. 
That bracket is easily deduced from the Poisson brackets of the components of 
L and e indicated above: for each energy level e, the map which associates to 



the Hamiltonian vector field with Hamiltonian 

restricted to the energy level e of the phase space, must be a Lie algebras homo- 
morphism. This remark leads to the formula 

[{ui,ut), (vt,V2)} e = y-Ul XVi + ^J^"2 X V2, -ui xv2+utxvt^J . 

By gluing together all the fibres tx^ over all points e G R, we get on R x 
£ x ~£^ a Lie algebroid structure, whose dual is the energy-mom entum space . The 
above formula gives the bracket of two smooth sections e I-* (iti(e), U2(e)) and 

e i y (y\{e),V2{ej) of that Lie algebroid. For each e G R, [(ui, u%), (vi, (e) 
only depends on the values taken by (1*1,1*2) and (vt,V2) at point e, because the 
considered Lie algebroid has a zero anchor map. 

3.11. Action of a Lie algebroid on a symplectic manifold. Let tta '■ A — > B 

be a Lie algebroid [3, 17] whose base is a smooth manifold B and whose anchor 
map is denoted by p : A —> TB. The bracket of two smooth sections s\ : B — > A 
and s 2 : B — > A of tta is denoted by {si, S2}. Let ttm ■ M — > B be a surjective 
submersion of a smooth manifold M onto the base B. An action of the Lie algebroid 
it a '■ A — > B on the fibered manifold ttm '■ M — > B is a map s 1— > X s which associates, 
to each smooth section s : B — > A of tta, a vector field X s on M, in such a way 
that for each smooth section s of tta and each smooth function / : B — > R, 

Xf s = (/ o tt m )X s , Tttm o X s = po SOTT A , 

and that for each pair (si,S2) of smooth sections of tta, 

X Sl +s 2 = + A S2 , [A Sl ,A S2 ] = A| S1 S2 } . 

Let us now assume that the manifold M is endowed with a symplectic form uj. 
The vector fields X s on M associated to smooth sections s of the Lie algebroid 
tta '■ A — > B cannot all be Hamiltonian, because if for some choice of the section 
s, X s is Hamiltonian, the vector field Xf s , associated to the section fs, where 
/ : B — > R is a smooth function, will not in general be Hamiltonian. So one may 
wonder what should be a reasonable definition of a Lie algebroid Hamiltonian action 
on a symplectic manifold. The answer is suggested by the example of the Kepler 
problem. 

In the last subsection, we described the Lie algebroid tta ■ A — >• B of infinitesimal 
symmetries of the Kepler problem: A = M. x£ x £ , B = R, tta '■ A — s- B is the 
projection on the first factor, the anchor map is the zero map A — > TB. The vector 
bundle tta '■ A -> B being trivial, each pair of vectors (mi,M2) G ^x? can be 
considered as a (constant) section ^ of our Lie algebroid. We defined the 
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action of our Lie algebroid on the phase space of the Kepler problem by taking, as 
vector field X S( _^ _^ associated to the section s^ ^, the Hamiltonian vector field 
with Hamiltonian 

This condition unambiguously defines our Lie algebroid action (and, simultaneously, 
the bracket composition law of its sections) since the module of its smooth sections 
is spanned by the constant sections. 

It would be interesting to see whether such a construction can be extended for 
more general Lie algebroids, with a nonzero anchor map. 



3.12. The S map. The main disadvantage of Moscr's rcgularization method is 
that it handles separately each energy level. This disadvantage may be partially 
removed by the following procedure, which allows to handle together all negative 
(resp., all positive) energy levels. However, negative, positive and zero energy levels 
still cannot be handled together with that procedure. 

Let us consider the quadric Q p , defined by 

h 2 + ((x 2 + y 2 + z 2 ) = (p 2 , 

where p may take any positive value, and let Qr be the quadric defined by 

h 2 + Q{x 2 + y 2 + z 2 ) = (R 2 , 

where R is a fixed positive quantity. To each point Mr £ Qr, we associate the 
point M p £ Q p such that 

P_ 

R 



OM p = £ OM R , 



We lift this diffcomorphism Qr h- Q p to the cotangent bundles, and we get a 
symplcctic diffcomorphism T p : T*Qr — > T*Q P . This symplcctic diffcomorphism 
associates to each pair (Mr, Wr) made by a point Mr £ Qr and a vector Wr 
tangent to Qr at that point, the pair (M p ,W p ) made by a point M p £ Q p and a 
vector W p tangent to Q p at that point: 

p —— * — > R 



OM p = % OM R , W p = - Wr . 
K p 

For each p > 0, we compose the symplectic diffeomorphism T p : T*Qr — > T*Q p with 
the symplectic diffeomorphism S p : T*(Q p \{N p }) -> T*£ built in Subsection 3.7. 
We obtain a family, indexed by p > 0, of symplectic diffeomorphisms S Pi r = S p oT p : 
T*(Qr\{Nr}) -> T*£: 



R-h, 



Op Rl with h 2 R + (\\Op R \\ 2 = R 2 



R - h R ^ , W hR 



r = W 3R H Op 



where we have set 



OM R = Op R + h R H , Wr = W 3R + W hR efi 
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The formulae for the inverse transformation are 



Our 



2Rp 

P 2 + (P : 



h R = R 



W 3R 



p 2 -Cp 2 
+ Cp 2 
<p 2 



V 



2Rp 



with p = = ^Jpl +p 2 y +pl 

-llp- f ^ 



hR 



R 



The diffeomorphism S R sends the subset of T*(£\{0}) on which the energy is 



2m 



into the subset of T*Qr on which 



km 2 
~Rp 



Therefore, if 



pi and pi arc two distinct possible values of p, the images by S~ R of the energy 

(p 2 Cp 2 

level EC 1 ? ,~$) = -, and by S~ 1 R of the energy level EC^,~0) = -, are 

2m P2, 2m 

disjoint. By restricting each map S~ R to the subset of T*(£\{0}) on which the 

Cp 2 

energy is E(r,p) = — — — , and by gluing together these restricted maps for all 

2m 

possible values of p, we obtain a unique diffeomorphism S^ 1 from the open subset 
of the cotangent bundle T* {£\{0}) minus the zero section on which the energy 
E(r,p) is negative if £ = 1, positive if £ = —1, onto T*(Qr\{Nr}) minus the 
zero section. This diffeomorphism, built by gluing together pieces of the symplectic 
diffeomorphisms S~ R for all values of p > 0, is no longer symplectic! It is given by 
the formulae, in which we no longer write the subscript i?, 



Op 



Ry/ 'C,r{2m 2 k — rp 2 ) ^ 

^ ___9 P ' 



h = 



R(rp 2 



m 2 k 
m 



2 k) 



W 3 = C 



w h = c 



m 2 k 
m 2 k~^ - 



?*(~r^. ~pf)~pf 



R v /(r(2m 2 k-rp 2 ) ' 
R 



We have set 



OAI = Op + het 
The transformed Hamiltonian (S~ 1 )*E = E o S is 

(k 2 m 3 



w h ei 



H = EoS 



2R 2 



the 



Ck 2 tn 3 

Up to the constant factor — — , it is the Delaunay Hamiltonian [4, 5], i.e.. 

Kepler Hamiltonian in Delaunay coordinates. This property is related to the fact 
that Delaunay variables are action-angle variables for the Kepler problem. 

In the rest of the paper we write Q for Q R . The north pole of Q is denoted by 
N instead of Nr. 

Let Xe be the Hamiltonian vector field on T*(£\{0}) with Hamiltonian E, i.e. 
the Hamiltonian vector field of the Kepler problem, and let Xh be the Hamiltonian 
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vector field, on T*(Q\{A}) minus the zero section, whose Hamiltonian H = E o S 
is given by the above formula. The direct image (S~ 1 )*(Xe) of the vector field Xe 
by the diffeomorphism 5 _1 is not equal to Xh, since S -1 is not symplectic. A short 
calculation leads to its expression, 

{S- 1 ) tf {X E ) = gX H , with 9 = j^- 

The function g is smooth on T*(Q\{JV}) minus the zero section, and becomes 
singular when h = R, i.e. on the fibre over the north pole N. We see that the 
map S^ 1 sends the Hamiltonian vector field of the Kepler problem to a conformally 
Hamiltonian vector field, with H as Hamiltonian and with g as conformal factor. 

3.13. The flow of Xh- Unlike the conformal factor g, which becomes singular on 
the cotangent space to Q at the north pole, the Hamiltonian H is smoothly defined 
on the whole cotangent bundle T*Q minus the zero section. The associated Hamil- 
tonian vector field Xh is complete. Its flow, defined on R x (T*Q\{zero section}), 

<S>x H (s, (OM{o),W0j)) = (OM(s),W^)) 

has slightly different expressions for £ = 1 (negative energy) and for £ = — 1 (positive 
energy). Observe that ||W|| is a first integral of Xjj- 

OM(s) = cos(As) OM(0) + = , sin(As) W(0) 
For C = l(^<0), - "^W" 



For C = -1 (E > 0) , 



Wc have set A 



W(s) = -ID^Jl s in(As) OM(o) + cos(As) W(0) 
R 

OM{s) = cosh(As) OM(o) - = > sinh(-As) W{o) , 

\\W(0)\\ 

W(s) = - JI^Mlj sinh(-As) OM(o) + cosh(As) W(0) 
R 



R 3 \\W{0 



3.14. A symplectic diffeomorphism. At the end of Subsection 3.12, we have 
seen that the image by S*" 1 of the Hamiltonian vector field Xe of the Kepler problem 
is a conformally Hamiltonian vector field gXn defined on an open dense subset 
of T*Q\{zero section}, whose Hamiltonian H is smoothly defined on the whole 
T*Q\{zero section}; the corresponding Hamiltonian vector field Xh is complete. 
Moreover the conformally Hamiltonian vector field gXn is also Hamiltonian, with 
Hamiltonian H, not for the canonical symplectic form of T*Q, but for the pull-back 
by S of the canonical s ymp l ectic form of T*£. Wc have seen in Subsection 3.4 that 
for every solution 1 1— > (r(t),p(t\) of the equations of motion of the Kepler problem, 




p(t).r(t)-2E(r(t),p(t))t \ = _J_ 
mk J r(t) 

The conformal factor g, whose expression is given at the end of Subsection 3.12, 
may also be expressed as 

mk / 1 
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Therefore 



p(*: 




i - 2E(r(f),p(t] 




2£(^,^ 


) 



All the integral curves of the conformally Hamiltonian vector field gXu can be 
written 

where i i-» (r(i) , p(t)) is a solution of the equations of motion of the Kepler problem. 
We see that the pull-back by (idu, S) of the function, defined on the product with 
K of the open subset of the phase space of the Kepler problem on which the energy 
E is negative if C = 1) positive if £ = — 1, 

has the properties of the function a of Theorem 2.4 of Subsection 2.3. That The- 
orem can therefore be applied. It proves that by composing the symplectic diffeo- 
morphism S^ 1 with the flow $>x H of the Hamiltonian vector field Xh, for suitably 
chosen values of the independent variable s, we get a symplectic diffeomorphism 
from the phase space of the Kepler problem, restricted either to negative, or to pos- 
itive energies, onto an open subset of T*Q. The explicit expression of this symplectic 
diffeomorphism is 



(f,^) 1 ^ (0AI(s),W(sty , with 



2E(~ 



The quantities (OM(s),W(s)) are given, as functions of (OM(0), W(0)) , by the 
formulae at the end of Subsection 3.13. We have to set 



OM(0) = Ofi + hU , W(0) = w 3 + w h r h . 

The formulae at the end of Subsection 3.12 give the expressions of O/i, h, W3 and 
Wh as functions of r and p. Finally we see that the symplectic diffeomorphism 
so obtained is that of Gyorgyi, Ligon and Schaaf. R. Cushman and L. Bates, in 
chapter II of the book [4], offer a detailed discussion of its properties, notably its 
behaviour with respect to the Liouville 1-forms of the cotangent spaces of the Kepler 
configuration manifold and of the 3-dimensional sphere (or hyperboloid) . 

4. Conclusion and perspectives. The properties of the diffeomorphism from 
the phase space of the Kepler problem to the cotangent bundle of a sphere or 
of an hyperboloid are explained, in Subsection 3.12, from a point of view other 
than that used by Cushman and Duistcrmaat [5]. Our explanation is founded on 
a very natural property of conformally Hamiltonian vector fields. However, the 
procedure we used, as well as those used by Gyorgyi [9], Ligon and Schaaf [15] 
remains unsatisfactory in the fact that it handles separately negative and positive 
energy levels. The space of motions of the Kepler problem is connected, since 
by varying slowly the energy level, an elliptic motion can be transformed into a 
parabolic, then into an hyperbolic motion. By completely different methods, J.- 
M. Souriau [21] proposed a global description of the manifold of motions of the 
Kepler problem, of its rcgularization and of its global and infinitesimal symmetries. 
It should be possible to do something similar at the level of the phase space instead 
of at the level of the manifold of motions. 
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The fact that the Kepler vector field is Hamiltonian with respect to a symplectic 
form and conformally Hamiltonian with respect to another syplectic form is re- 
markable; it is probably related to the complete intcgrability of that vector field, 
as suggested by an anonymous referee. In [16], A. Maciejewski, M. Prybylska and 
A. Tsiganov have used conformally Hamiltonian vector fields within the theory of 
bi-Hamiltonian systems, to build completely intcgrablc systems. 

Wc have shown that the infinitesimal symmetries of the Kepler problem form a 
Lie algebroid rather than a Lie algebra. Actions of Lie algebroids on symplectic 
manifolds is a subject which seems to us interesting; it should be nice to have other 
examples of Hamiltonian vector fields whose infinitesimal symmetries form a Lie 
algebroid, maybe with a non-zero anchor map. 

A. Douady and M. Lazard [6] have shown that Lie algebroids with a zero an- 
chor map are intcgrable into Lie groupoids; it would be interesting to write down 
explicitly the action of a Lie groupoid with the groups SO(4), SE(4) and SO(3, 1) 
as isotropy groups on a regularized phase space of the Kepler problem. 
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